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O ■ Abstract 

Sterile neutrinos with a mass around the keV scale are an attractive particle physics 
candidate for Warm Dark Matter. Although many frameworks have been presented 
in which these neutrinos can fulfill all phenomenological constraints, there are hardly 
any models known that can explain such a peculiar mass pattern, one sterile neutrino 
at the keV scale and the other two considerably heavier, while at the same time being 
compatible with low-energy neutrino data. In this paper, we present models based 
on the Froggatt-Nielsen mechanism, which can give such an explanation. We explain 
how to assign Froggatt-Nielsen charges in a successful way, and we give a detailed 
discussion of all conditions to be fulfilled. It turns out that the typical arbitrariness 
of the charge assignments is greatly reduced when trying to carefully account for 
all constraints. We furthermore present analytical calculations of a few simplified 
models, while quasi-perfect models are found numerically. 
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1 Introduction 



One of the most intriguing problems in today's particle physics and astrophysics is the 
identity of the mysterious Dark Matter [1] observed in the Universe [2]. It is known since 
a few years that computer simulations of small scale structure formation are in favour of 
so-called "warm" Dark Matter (WDM), with a mass of roughly 1-2 keV [3j H], which is an 
intermediate case between the standard paradigm of "cold" (non-relativistic) and "hot" 
(relativistic) Dark Matter, the latter being strongly disfavored by structure formation. 
Furthermore, there are model-independent data analyses which also seem to point to the 
keV scale [5]. 

A very interesting framework providing a WDM candidate motivated by particle physics 
is the so-called z/MSM [6], which extends the Standard Model (SM) of particle physics by 
three right-handed (sterile) neutrinos, one of which has a mass at the keV scale, whereas 
the other two are considerably heavier. However, it turns out that, in the standard thermal 
production, the Dark Matter would be overproduced, which makes it necessary to rely on 
non-thermal production instead [TJ EJ [9]. An alternative is provided by embedding the 
keV sterile neutrino in a gauge extension of the SM and correcting the abundance by 
allowing for sufficient entropy production in the decay of the two heavier sterile neutrinos, 
which has been exemplified in a Left-Right symmetric framework [TO]. Note, however, 
that this last proposal required a seesaw type II situation in order not to be in conflict 
with the experimental and observational constraints. In addition, keV sterile neutrinos can 
also appear in the frameworks of the scotogenic/inert Higgs doublet model [HI [12] or of 
composite neutrinos |13j . 

The common feature of all these proposals is, however, that they do not yield an expla- 
nation of the required mass pattern of sterile neutrinos: They can be seen as very useful 
frameworks providing all tools to make concrete predictions, and the corresponding pa- 
rameters can assume values leading to full agreement with data. Nevertheless, in the view 
of model building, they only assume the correct mass pattern to be present. 

Up to now, to our knowledge, only two classes models have existed which could yield an 
explanation of the sterile neutrino mass pattern. The key point is to achieve a strong mass 
splitting in the right-handed neutrino sector, as only the lightest sterile particle should 
have a mass in the keV range, while the other two must be considerably heavier [TO] . This 
leads to different schemes for shifting a certain initial mass spectrum, two of which are 
depicted in Fig. [TJ The first class of models is based on a flavour symmetry which forces 
one sterile neutrino to be strictly massless. This idea has first been discussed in Ref. [TT] 
for a specifically defined lepton number symmetry, and has recently been investigated in 
the context of a L e — L M — L T flavour symmetry [15] . Both symmetries can explain such 
a peculiar pattern of right-handed neutrino masses easily: The trick is to make use of 
the different scales of symmetry breaking and symmetry preserving terms to generate the 
necessary hierarchy. Soft breaking of the symmetry will lift up the mass of one sterile 
neutrino, which would have been strictly massless in the limit of restored symmetry, to 
the keV scale, see left panel of Fig. [TJ for the L e — — L T case. The breaking will 
furthermore break up the predicted exact mass degeneracy in the light neutrino sector, 
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Figure 1: The mass shifting schemes of the two classes of models that had already existed: 
The left panel displays the shifting due to soft breaking of a global L e — — L T flavour 
symmetry, while the right one shows the exponential suppression in a Randall- Sundrum 
framework. 



which would otherwise contradict the data. The second real model [16] exploited the 
exponential factor in a Randall-Sundrum [T7] like framework to obtain a large mass splitting 
from very moderately tuned parameters, cf. right panel of Fig. [TJ This model has the 
advantages of naturally explaining the existence of one sterile neutrino at the keV scale, 
while the other two could have masses of around 10 11 GeV or even heavier, and of at the 
same time insuring that the seesaw mechanism JT8J [19j EH IZEJ [22] works, in spite of the 
existence of a relatively light sterile neutrino. On the other hand, this model involves the 
assumption of a suitable UV-brane, which is intrinsically hard to probe. 
In this paper, we want to present a third alternative to explain the existence of a keV 
scale sterile neutrino, by making use of the Froggatt-Nielsen (FN) mechanism [23]. This 
mechanism is well-known to be capable of explaining very strong mass hierarchies, e.g. in 
the quark sector [211 [2S1 [26], which is just what is needed for having one sterile neutrino 
at the keV scale whereas the other two are heavier by at least six orders of magnitude. 
Our goal is to find the minimal assignments needed to explain the pattern in the sterile 
neutrino sector, while being in full agreement with the rest of the lepton data. The key 
point of the FN mechanism is that at most one fermion generation (namely the third one, 
if at all) obtains a mass by the standard Higgs mechanism, or by the existence of a direct 
mass term for gauge singlets. The other generations, however, will only receive higher- 
order contributions from multiple seesaw-like diagrams, which leads to the typical cascade 
or waterfall structures |27] . This is enforced by having an additional Abelian [7(1)fn 
flavour symmetry under which the first, second, and third generations have charges with 
decreasing absolute values. Furthermore, such a framework requires a heavy sector of new 
fermions that is usually not specified in details, as well as so-called flavon fields, which 
are SM singlets but charged under the £7(1)fn- These fields obtain vacuum expectation 
values (VEVs) in order to break the symmetry in a phenomenologically suitable way [28] . 
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Integrating out the heavy fields finally leads to suppression factors by some power of a 
small parameter A, which is usually taken to be of the order of the Cabibbo angle [29] . 
Note that, shortly before this work was completed, Ref. [30] appeared, where a Froggatt- 
Nielsen mechanism, together with a Z3 symmetry, is used in the context of an A4 flavour 
symmetry to generate a hierarchy for the charged lepton masses and to regulate the mass of 
a sterile neutrino to be around the eV-scale. The authors also comment on the possibility 
of exploiting this model for an explanation of a sterile neutrino with a keV mass, however, 
without investigating this route explicitly. In their case, however, the Froggatt-Nielsen 
mechanism is merely a small addendum to the existing A4 model, whereas we consider the 
mechanism as starting point to develop fully working models. 

The structure of this paper is the following: In Sec. |2] the main features of the Froggatt- 
Nielsen mechanism are summarized The possible choices for the FN charges of the right- 
handed neutrinos are analyzed in Sec. [3j After that, in Sec. HJ we give a detailed discussion 
about all requirements to be fulfilled, where we will also explain why certain frameworks 
are incompatible with the FN mechanism when aiming at a description of keV mass sterile 
neutrinos. In addition, we will also discuss why certain constraints that might be prob- 
lematic for a higher right-handed neutrino mass scale are practically of no relevance to our 
case. Finally, in Sec. |5l we present analytical calculations of some promising scenarios, as 
well as fully working numerical models, before concluding in Sec. EJ Explicit expressions 
for the diagonalization matrices used can be found in Appendix A. 

2 The Froggatt-Nielsen mechanism 

The Froggatt-Nielsen (FN) mechanism [23] is maybe one of the best possibilities to explain 
strong hierarchies between fermion masses. Furthermore, there exist successful explicit 
applications of this mechanism to the neutrino sector, see, e.g., Refs. [261 [31], an d Refs. j32j 
1331 [3U [35] for a systematic parameter space scan. 
We would like to use the FN mechanism to explain a hierarchical spectrum (HS) in 
the sterile neutrino sector of the following type: Mi ~ O(keV) and M 2 ,M 3 > O(GeV). 
Denoting the FN flavon field by 0, and the cut-off scale at which the heavy sector of the 
theory is integrated out by A, we can write the Lagrangian as: 

/Q\h+fj _ /Q\9i+fj 

Aeptons = -Y: j e^HL jL I -J + h.c.-Y%N iR HL jL (-\ + h.c. (1) 

-- N iR Ml (N jR f (-] + h.c. - - Y? (L lL )c (*a 2 A) L jL (-) + h.c. , 

where H is the Standard Model Higgs and H = ia2H* is its charge conjugate (see Tab. [1] 
for the definition of all the other fields). The FN flavon field can acquire a VEV (G), with 
A = being a small quantity of the order of the Cabibbo angle: A ~ 0.22 |29j . 
Note that, in principle, FN charges could be positive or negative, which would correspond 
to integrating out heavy particles or anti-particles, respectively. This would again lead to 
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Table 1: SU{2) and family charges, with i = 1,2,3. 



suppression factors due to having a higher mass scale involved. Accordingly, the decisive 
quantity for the suppression is actually the absolute value of the sum of FN charges. After 
the flavon field acquires a VEV and after electroweak symmetry breaking, the charged 
leptons and the Dirac neutrino mass matrices are given by 

(yll A |fci+/i| yl2^|fci+/ 2 | yl3 A |fci+/ 3 |\ 
yJi^+M y^2 A |fc 2 +/ 2 | y^3 A |fc 2+ / 3 | ; ( 2 ) 
yjU A |fc 3 +/i| yj32 A |fc 3 +/ 2 | yjS3 A |fc 3 +/ 3 | / 

(y^i \\ai+fi\ y^ 2 A l9l+/21 Y^ 3 A |91+/3| \ 
y|i A | S2 +/i| y|2 A | 92 +/ 2 | y|3 A | 52 +/ 3 | ? ( 3 ) 
yJi A | S3 +/i| yJ2 A | 93 +/ 2 | yJ3 A | 93 +/3iy 

where t> is the VEV of the Higgs doublet and fi, ki, and gi are the FN charges of the 
left-handed lepton doublets, the right-handed charged, and the right-handed neutrinos, 
respectively, see Tab. [TJ 

Now let us have a look at the right-handed neutrino sector. Since the corresponding mass 
matrix Mr has to be symmetric, it will have the following structure: 

/M^X^ M^X^+92\ Mj f 3 A l3i+ 9 3|\ 
Mr = I • MfA^I MfAl^l . (4) 
\ • • MfAl^l / 

Applying the type I seesaw formula then leads to the following structure: 

'aiAl 2 *l h\^ + M dAl^+^l 
m[ = -m T D M^ 1 m D = \ • dxX^ a X^ 2+f ^ \ , (5) 

where the parameters a±, b±, c\, d±, e±, fx depend on the parameters present in the mass 
matrices mo arid Mr. In the case of a type II seesaw scenario, one needs the left-handed 
Majorana neutrino mass in addition, 

/yn A | 2/l | y L 12 A |/ 1+ / 2 | yl3 X \f 1+ f 3 \\ 

m L = v A \ • FfAl^l Y^XM+M }, (6) 
V • • FfAl^l / 

with va being the VEV of the triplet A. The type II seesaw neutrino mass matrix will be 
given by 

'a 2 Al 2/l l 6 2 Al /l+/2 l c 2 Al /l+/3 l 



m 



m L - m T D M^ 1 m D = m L + m[ = ( • d 2 Al 2/2 l e 2 Al /2+/3 l | , (7) 

f 2 X™ 
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where the parameters 02, 62, C2, d^, e<i-> ji depend on the parameters present in the mass 
matrices m L , m D , and M R . 

Note that Eqs. (JSJ) and (0) proof that the seesaw mechanism is, in this framework, not 
spoiled by the presence of a keV neutrino: Any global U(l) charge would cancel out in the 
seesaw formula, and so do the C/(1)fn charges of the right-handed fermions. These charges 
are, however, the only connection to the keV scale, since they lower the natural scale M Q 
of the right-handed neutrino mass by introducing certain powers of A. Accordingly, in the 
seesaw formula, one is always guaranteed to divide by a relatively large value M , thereby 
saving the seesaw mechanism. In fact, since the FN charges of the left-handed lepton 
doublets are still present in the seesaw formula, the suppression mechanism for the light 
neutrino masses might even be amplified. 

Furthermore one can see from Eqs. (jSJ) and (JTj) that the structure of the neutrino mass 
matrices is precisely the same for type I and type II seesaw scenarios, which is another 
consequence of the right-handed neutrino charges canceling out. For this reason, we will 
not consider type II seesaw scenarios any further in this paper, since the results for such 
scenarios can be trivially recovered from the seesaw type I results, by the simple transfor- 
mations a\ — > a,2, bi — > 62, and so on. We will, however, mention differences or additional 
features that would appear in the type II seesaw case, where appropriate. 

Before starting to apply the FN mechanism to explain keV sterile neutrinos, we also want 
to stress a potential problem of this method: The Froggatt-Nielsen mechanism intrinsically 
involves a high energy sector that is not specified further. This is a problem not only of 
this particular mechanism, but of practically all flavour symmetries that must be broken 
in a phenomenologically acceptable way. Such a breaking always involves a scalar (flavon) 
sector, which is assumed to have suitable properties. In addition, the FN mechanism 
also involves heavy fermions, which are integrated out to lead to the suppression factors 
in Eq. (jTJ) . Such high energy sectors, although of no practical relevance for low energy 
experiments, could potentially become important for very high energies, and hence in 
particular in the early Universe. This problem, however, is strongly mo del- dependent, and 
far beyond the scope of a conceptual paper like the one presented here. 

We nevertheless want to point out that the additional interactions generated by this high 
energy sector could, e.g., cause the keV sterile neutrinos to be in thermal equilibrium at 
early times. Although certain production mechanisms of keV sterile neutrino DM require 
the sterile neutrinos to never enter thermal equilibrium (see, e.g., Refs. [6j [7J [8j [9]), there 
is also the alternative possibility to first produce them thermally and then dilute their 
abundance by sufficient entropy production [TO]. The lesson to learn is that one should 
be careful when applying the FN mechanism to certain scenarios: We will discuss many 
accompanying problems from the particle physics side later on in Sec. HI However, one has 
to keep in mind that also certain astrophysical scenarios could lead to further restrictions 
or incompatibilities. 
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3 FN charges for the right-handed neutrinos 



The first question we want to answer is how to assign FN charges to the right-handed 
neutrino fields. The most important constraint to keep in mind is the strong hierarchy 
in the right-handed (sterile) neutrino sector: In a framework as in Ref. [10], where the 
lightest heavy neutrino is supposed to have a mass M x of a few keV, the second to lightest 
(heavy) neutrino must at least have a mass M 2 of order GeV, due to the requirement of 
sufficient entropy production. This condition does, however, not constrain the largest mass 
M 3 , which can hence be similar to or even much larger than M 2 . In any case, we require 
a hierarchy of at least six orders of magnitude between M\ and M 2 . 
If we suppose that all the coefficients in the matrix M# are of the same order, then the 
explanation of the HS in the sterile neutrino sector should come from the FN charges. 
Considering the FN parameter A to be of the size of the Cabibbo angle [29] (A ~ 0.22) and 
9i > 92 > 93, the minimal conditions to be fulfilled are: 

f 9i > fi'ilmin with #i| min = 5-2 + 3, 
\ 92 > 92 1 mm with g 2 \ min = g 3 . 

In the minimal case, g 2 = g 3 = gi — 3, we would obtain a mass spectrum of the following 
type: M l ~ C(keV) and M 2 ,M 3 ~ C(GeV). For g 2 = g 1 - 3 and g 3 < gi, the mass 
spectrum would instead be given by: Mi ~ C(keV), M 2 ~ C(GeV), and M 3 > C(GeV). 
Finally, if g 2 < g\ — 3 and g 3 < gi, both Mi and M 3 have a mass greater than 0(GeV). 
Since we want to stick to minimal choices of FN charge assignments, a reasonable condition 
would be to set g 3 = and g± = gi + 3. Accordingly, we choose two minimal scenarios that 
we will investigate further: 

• Scenario A: (g u g 2 , g 3 ) = (3,0,0), 

• Scenario B: (g u g 2 , g 3 ) = (4,1,0). 

These scenarios lead to mass eigenvalues M 12 3 which obey the hierarchies Mi ~ 10~ 6 M 2 3 
and Mi w 10 _6 M 2 ~ 10 _8 M 3 , respectively. In Tabs. [2] and [31 we show some examples of 
allowed textures that lead to the desired HS, considering two, three, or four independent 
parameters in the mass matrix Mr. 

As explained before, one could increase the splitting of the mass eigenvalues by assuming 
stronger hierarchies in the FN charges. Furthermore, one could in principle also assign a 
non-zero g 3 , which would decrease the values of the masses compared to some characteristic 
scale that could, e.g., be generated by the VEV of some scalar field. 

In Fig. [21 we have schematically depicted the general effect of the FN mechanism: A certain 
mass scale Mq is multiplied by powers of A that depend on the fermion generation. These 
factors lead to suppressions of the physical mass eigenvalues. In general, FN assignments 
are very well suited to explain strong hierarchies, which is why we ultimately chose to 
investigate this framework. 

In the following, we will investigate only the two exemplifying Scenarios A and B, and we 
will show how to implement their assignments in a more complete model. 
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Table 2: Examples of Mr textures for the Scenario A that lead to a HS. We have assumed 
all the parameters to be greater than zero. 
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Table 3: Examples of Mr textures for the Scenario B that lead to a HS. We have assumed 
all the parameters to be greater than zero. 

4 Requirements, No-Go's, and amenities 

In this section, we will explain in a careful way why several FN scenarios are not able to 
successfully produce a neutrino sector compatible with the data, or at least have problems 
with it. In fact, even though the FN framework seems to involve quite some freedom, 
the way how the charges are combined is nevertheless quite restrictive: By fixing only 9 
charges (lepton doublets, right-handed electrons, and right-handed neutrinos), we predict 
the magnitudes of (9 + 9 + 6) = 24 (type I seesaw) or (9 + 9 + 6 + 6) = 30 (type II seesaw) 
entries in the matrices M e , mo, Mr, and rriL, which would otherwise be independent. 
Accordingly, there can be many allowed theoretical and phenomeno logical scenarios that, 
when implemented within a FN framework, may lead to contradictions. In the following, we 
will discuss several situations which can lead to problems, as well as certain requirements 
that have to be fulfilled, thereby reducing step by step the arbitrariness involved in the 
FN models. Furthermore, we will point out why certain constraints that are typically 
problematic for FN inspired models do not apply in our case. As it will turn out, although 
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Figure 2: The mass shifting scheme of Froggatt-Nielsen models. 



FN charge assignments might seem relatively random at the first sight, they are in fact 
restrictive enough to disagree with several frameworks, and in particular they disagree 
with the Left-Right symmetric framework for keV sterile neutrino Dark Matter proposed 
in Ref . [TU] , as well as with the bimaximal mixing from the neutrino side proposed in this 
context in Ref. [T3|. Turning this round, it will be easy to rule out the FN framework if 
evidence for one of the contradicting scenarios is found. 

Note that one could, in principle, apply additional flavour symmetries to force some entries 
in the mass matrices to zero. This method could alter all conclusions in this section, but 
such an approach is, on the other hand, far away from being minimalistic. In particular, 
the FN mechanism might not even be necessary anymore in such cases if, e.g., some VEV 
hierarchies are imposed or if the symmetry is softly broken. Since we want to stick to FN 
situations, however, we do not consider such extended scenarios in this paper. 

4.1 No bimaximal neutrino mixing 

Let us now investigate the compatibility of certain scenarios with the FN framework. One 
easy scenario that leads to a tri-bimaximal PMNS matrix [36] is given by a bimaximal 
neutrino mixing and an opportune form of the charged lepton mixing, as discussed in 
Ref. [37]. When trying to implement this scenario with a FN mechanism, one however 
encounters problems. Indeed, the assumption of bimaximal neutrino mixing poses strong 
constraints on the form of the light neutrino mass matrix, which can, most generall}0, look 

1 Actually, there could be a second contribution that is exactly proportional to a unit matrix. Such a 
contribution, however, could never be explained by the FN mechanism, which can only explain orders of 
magnitude, but no exact equalities. Accordingly, this even more general case is not of practical relevance 
for FN inspired models. 
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like [38]: 

/O A B' 

M bi-max = . 

V- • o, 

Comparing this form with Eqs. (jSJ) and (J7J), we can see that the FN charges (fx, f 2 , fa) of 
the left-handed leptons need to be such that |/i + /2,3| is small, while all other combinations, 
and in particular |2/j|, have to be large. This enforces large and more or less equal absolute 
values of all fa's, while simultaneously sign(/i) = — sign(/ 2 ,3) must be fulfilled. This is a 
very specific choice that cannot be brought in agreement with any of the possible forms 
of the charged lepton mass matrix M e , listed in Ref. [38J, which could lead to a leptonic 
mixing in agreement with the experimental values. We have checked that these conclusions 
are not significantly altered by the use of two flavon fields. 



4.2 No Left-Right Symmetry 

As we have just seen, it is relatively easy to impose requirements on the FN charges that are 
too restrictive to be fulfilled. This problem originates from the fact that it is non-trivial to 
obtain large mixings, as required for leptons, from a (non-lopsided) FN framework [39]. As 
already seen in Sec. [2j the symmetric form of the Majorana mass matrix essentially leads to 
some cascade-like structure at best, which can nevertheless lead to tri-bimaximal leptonic 
mixing in case that also the charged lepton mass matrix has a cascade-like form [27]. It is 
exactly this last requirement, however, that is spoiled by Left-Right (LR-) symmetry. 
Imposing Li?-symmetry and including Higgs triplets [10] in order to accommodate for the 
type II seesaw situation that is required in that context [10] . the most general leptonic 
mass matrices must have the forms 

(m D )ij = vifij + v 2 gij, 

(M e )ij = vxgij + v 2 fij, 

{rriL)ij = V2v L hij, and 

{M R )ij = V2v R hij, (9) 

where v\ )2 are the Higgs doublet and vl,r are the Higgs triplet VEVs. However, in such 
a model the right-handed charged leptons e^R are grouped into doublets $ R of SU(2) R 
together with the right-handed neutrinos N iR [and hence their FN charges must be equal, 
(ki, k 2 , k 3 ) = (gi, g 2 , #3)], and the discrete Li?-symmetry dictates the equality of the abso- 
lute values of the FN charges between the left- and right-handed doublets, Q(^ l L ) = Q(^ 1 r) 
[and hence (fi, f 2 , f 3 ) = (gi,g 2 ,g3), if we consider only positive FN charges] [2B]. From 
these conditions, our two Scenarios A and B already fix the complete structure of the mass 
matrices. Hence, the most general forms for the charged lepton and for the light neutrino 
mass matrices in Scenarios (A,B) are given by 

/l A 3 A 3 > 4 \ /A 6,8 A 3,5 A 3 < 4 \ 

M e oc A 3 1 A ' 1 and m L oc • A 0,2 A ' 1 . (10) 
\A 3 < 4 A ' 1 1 / \. . 1 / 
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Since M e is a Dirac-type mass matrix, and since the above form is dictated directly by the 
choice of (g±, g^-, gal, we cannot avoid a large 11-element, which clearly spoils the demanded 
cascade structured In particular, none of the assignments that we will use later in the 
working examples fulfills the condition Q(^ l L ) = Q(ty l R ). 

4.3 More than one FN field 

It had been shown in Ref. [H] that, in a seesaw framework, it is difficult to obtain a large 
mixing angle scenario for leptons with only a single £7(1)fn, unless one relies on pseudo- 
Dirac scenarios, which essentially involves setting some elements of the mass matrices equal 
to zero |12] . 

A very interesting comparison between one or instead two flavon fields has been performed 
in Ref. [13] for the quark sector in an SU(5) Grand Unified Theory (GUT) inspired scenario, 
and it has been extended in Ref. [26] to the lepton sector. These references agree that, 
in order to have models which can be treated analytically, only real entries in the mass 
matrices should be investigated. This, however, will not allow for CP violation, since all 
phases could be trivially rotated away. A way to accommodate non-trivial CP phases is 
to introduce two flavon fields 0^2 rather than only one, and to require the VEVs of these 
two fields to have a relative phase. In addition, 0i2 must have opposite charges under an 
auxiliary Z2 parity in order for this phase to survive. A further problem of models with 
only one FN field is that they normally lead to small atmospheric neutrino mixing [26], 
and are thus incompatible with the data. For these reasons, we are going to analyze in 
Sec. [5] the case of two flavon fields, instead of having only one. 

4.4 Anomaly cancellation in SU(5) 

The FN charge assignments that we are going to use are coming from an SU(5) GUT 
scenario [H], as mentioned in Ref. [26] where, however, no detailed explanation for this 
form was given. We will do this by having a closer look at the conditions that have to be 
fulfilled in order to guarantee the cancellation of dangerous anomalies. 
The charge assignments used by us can be easily understood by looking at the conditions 
that have to be fulfilled in order to guarantee the cancellation of anomalies. Such con- 
siderations are used, for example, in the models from Refs. [15| H6"] . Here we will follow 
the procedure outlined in Refs. [311117]. In an SU(5) GUT model, the right-handed elec- 
tron is situated together with the quark doublets and the right-handed up-like quarks in a 
10-representation, whereas the lepton doublet is grouped together with the right-handed 
down- type quarks in a 5-representation. We use the parametrizations of Ref. [31], where 
the authors define the FN charges of the quark doublets, of the right-handed up quarks, 

2 Our statement, however, only holds with certainty in a FN framework, and in particular for Majorana 
neutrinos, i.e., the light neutrino mass matrix is demanded to have a symmetric structure. Note that in 
the framework of Ref. [57] , the authors suggest Li?-symmetry to lead simultaneously to cascade structures 
in M e and m„, which can be fulfilled in a non-FN context. In the cases discussed here, however, this is 
not possible. 
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and right-handed electrons, respectively, as 

3 3 3 

qj = x + u , ^^Ui = x + 2u, ^^ei = x, (11) 

i=l i=l i=l 

while the lepton doublets and right-handed down quarks are 

3 3 

^2i i = y , J2 d i = y + v - ( 12 ) 

i=l i=l 

Anticipating the Assignments (1,2) to be used later on in Sec.0 cf. Tab.HJ we can determine 
the parameters x, y, u, v from the conditions 

k\ + k 2 + k 3 = 5 = x + u = x + 2u = x, 
/i + / 2 + / 3 = (l,4) = y = y + v, (13) 

which immediately lead to u = v = (as characteristic for assignments consistent with 
SU(5) [31]), and hence also to x = 5 and y = (1, 4). Since we will consider a case in which 
the SM-like Higgs (or other Higgses) do not carry any FN charge, we also need to fulfill 
= — z = z + u + v [31], which is no problem if z — 0. Then, we can easily satisfy the 
condition for anomaly cancellation: 

A 3 = A 2 = ~A X = ^[3x + Au + y + v] = -[3^ + k 2 + k 3 ) + h + f 2 + /a] = (8, 9.5), (14) 
while the A[ vanishes for both assignments, as demanded: 

A[ = + ^ -<% + l*-e*) = 0. (15) 

i=l 

Note that the FN charges {gi,g 2 ,gs) of the right-handed neutrinos do not appear in the 
conditions displayed in Eqs. ( IT3|) and ( IT4|) . They are, indeed, total singlets not only under 
the SM gauge group, but also under SU(5), and hence they cannot contribute to any gauge 
anomaly. This is the key point to be able to freely implement our right-handed neutrino 
Scenarios A and B. 

4.5 Difficulties with SO(10)? 

We have decided to analyze in detail an SU(5) inspired model. The reason for this is 
not only that the right-handed FN charges {gi,g 2 ,g-3) drop out of the light neutrino mass 
matrix, as they would in any type I or II seesaw model, but also that they are essentially 
unconstrained, since the right-handed neutrino iVj is a singlet l_j for each generation i in 
SU(5) @3|. In an 50(10) GUT [H], instead, the right-handed neutrino A 7 * is part of the 16^ 
representation together with all quarks and leptons of generation i [48]. This requirement 
would constrain the right-handed neutrino FN charges strongly, and it is therefore not 
clear if it is possible to find a realistic setup that reproduces all data correctly, while at 
the same time keeping a strong mass splitting in the right-handed neutrino sector. 
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4.6 Problems with democratic Yukawa matrices 



From a FN model, we generally expect a waterfall structure in the charged lepton and 
in the neutrino mass matrices [27]. Let us suppose that a specific FN model leads to the 
following matrices, which have just the form that we will also obtain for both our scenarios, 
cf. Sec. apart from overall factors: 



A 4 


A 3 








A 




A ? > 


A 2 




, m v oc 




1 




A 2 


A 






K A 








M\M e oc A 3 A 2 A , m u oc A 10. (16) 
V A 2 A 1 / \ X 1/ 

Then, these matrices will be diagonalized respectively by U e and by U u ~ U\Uh' 

/ 1 \ / 1 A A \ / 1 \ 

U e ~ 1 A I , ?7x — ( —A 1 I , Uh - 1 ] , (17) 




with U\ being the matrix that diagonalizes m„, and Uh the matrix that corrects the 
neutrino eigenvalues for inverted ordering. Note that the correction Uh is required due to 
the presence of two large and only one small mass eigenvalues, which cannot be realized 
in normal hierarchy. In this case the PMNS matrix, 




has a form which is not consistent with the data [19], as can be easily seen by calculating 
the value of # 13 . The key points of this result are a waterfall structure for the mass matrices 
and the presence of the matrix Uh, and thus of an inverted ordering scenario in the neutrino 
sector. 

This simple argument shows that in general we need to go beyond the democratic Yukawa 
coupling assumption in a FN scenario. There is an easy way to overcome this problem, 
namely by considering a slightly non-democratic Yukawa hypothesis. We will demonstrate 
this explicitly in Sec. I5.1.2l by considering specific SU (5) inspired models, while in Sec. 15.1.11 
we will analyze the democratic case and find indeed a PMNS matrix of the form of Eq. (|T8|) . 



4.7 No need for RGE running 

The mass matrices we obtain are actually only correct at a high energy scale, like the 
GUT scale, where the FN charges are imposed. Since, however, the neutrino observables 
are measured at a low energy scale, we have to evolve the neutrino masses and mixing 
parameters down to that scale by renormalization group equations (RGEs) [50l |5T| [52] , 
an effect that is often dubbed as "running". Sometimes, this step is not applied, and it 
is instead argued that this should be the reason why a certain model does not fit to the 
data [26], although this must not necessarily hold true. But in general the running has 
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to be taken into account, as there is no way to avoid it. We will, however, show in the 
following that running effects, although present, are fully negligible in our case. 
First, note that the running of the charged lepton Yukawa coupling matrix Y e and of the 
light neutrino mass matrix m v in a type I seesaw model are given by [50] . 



dY 

16tt 2 — - = Y e {D e Y^Y e + D u Y ] d Yd)+ (flavour diagonal), (19) 

(JjL 

16tt 2 ^ = (C e Y}Y e + C u Y ] D Y D ) T m y + m v (C e Y ] e Y e + C U Y ] D Y D ) + (flavour diagonal), 

where t = ln(/x//xo), fJ> is the renormalization scale, and /io is the reference scale (e.g. the 
GUT scale) at which the input information is imposed. The coefficients C eji , and D e v are 
numbers of 0(1). The flavour diagonal terms are not displayed, since they could only 
lead to an overall rescaling which is not important for mass ratios, but they will not affect 
the mixing. Note that we have neglected subtleties like threshold effects as they are not 
relevant to our argumentation. 

Now, since we are considering a very low-scale seesaw framework, where the right-handed 
neutrino mass scale could be as low as a few GeV, the seesaw formula m u = —m^M^rriD 
together with the definition of the Dirac neutrino mass, = yov, implies that the order 
of the Dirac Yukawa coupling must be as small as, e.g., yo ~ 10~ 5 for m u = 1 eV and 
Mr ~ 10 GeV. The dominant entry in Y e , however, is y T ~ 0.01. Accordingly one can, in 
Eq. (fT9|) . completely neglect Yd, and safely assume that the largest number on the right- 
hand side is about 0.01 2 = 0.0001. Dividing this number by the loop factor 167r 2 decreases 
it to a value that is even smaller than the electron Yukawa coupling. Accordingly, any 
flavour non-diagonal term on the right-hand sides of Eq. (1T91 will be small. Hence, the 
only effect the running can have is an overall scaling. Furthermore, due to the smallness 
of Yd, one can also neglect the running of the right-handed neutrino mass matrix [50J. 

We can show more explicitly that the mixing angles and phases do not undergo a con- 
siderable running. The correction to the mixing angles can be estimated to be at most 
about [50] 

A6~-^y^ln(^), (20) 

where £ is an enhancement or suppression factor that can be at most as large as £ ~ ^4 ~ 

25. Hence, from Eq. (|20|) . one can estimate the maximal correction to a mixing angle to 
be about 0.001, which is perfectly negligible. Similarly, the evolution of the phases is, for 
a general phase 0, roughly given by [5T] 

2 i 



A0~— ^ 2 -Mn ^ , (21) 

16tt 2 Qj V /' / Qj 

where (a = mi ~ mj is a function of the light neutrino masses. Obviously, this correction to 
the phase is also completely negligible, unless strong degeneracies in the neutrino masses 
lead to a very small Qj. Such a situation, however, is practically impossible to achieve 
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in FN models, since any sensible FN charge assignment will always introduce hierarchies 
rather than degeneracies, so that we do not have to consider the running of any phases. 
Accordingly, the only effect that running could have is an overall scaling of the mass 
matrices. Such scalings are, however, implicitly included in the prefactors of our mass 
matrices and, in particular, they will cancel out in mass ratios or ratios of mass squares. 
Furthermore, the running will practically not affect any mixing angles or CP phases. We 
have verified numerically for several examples that this is indeed the case. In fact, for 
Yukawa couplings that are not larger than about y T , numerical computations show no sign 
of running in the mixing angles in the absence of extreme degeneracies, even if we run over 
several orders of magnitude in energy |52j . 

4.8 Potential constraints from lepton flavour violation 

There is also a (seemingly unrelated) problem we want to comment on, which is, from 
a theoretical point of view, not necessarily connected to our models, but which might 
nevertheless show up in practice. This is the generic phenomenon that theories beyond the 
SM do not easily respect flavour, and will hence tend to lead to lepton flavour violating 
(LFV) reactions [531 ED- 

This problem arises in our case because we rely on a GUT framework which, in turn, 
requires a unification of the SM gauge couplings. However, with only the SM particle 
content, this unification does not happen [55]. In principle, one would not necessarily 
have to care about this problem in a FN framework, since the FN mechanism intrinsically 
involves the existence of a high energy sector that is not specified further, and which is 
assumed to have just the right properties as to make the FN mechanism work [23J. The 
easy way to go would be to simply assume this high energy sector, which is present anyway, 
to also be responsible for gauge coupling unification. 

However, in practice, one would like to have a high energy sector at hand that is specified, 
in order to make concrete predictions. One of the known ways to achieve gauge coupling 
unification is to assume the presence of supersymmetry (SUSY) [55], which is often con- 
sidered when talking about GUTs. But the introduction of SUSY also leads to problems 
connected to this theory, the prime example being generically large LFV effects [56] . 

These effects have been studied in the context of GUT-inspired FN models, like ours, in 
Ref. [39]. Even if a flavour diagonal universal slepton mass mo is assumed at the GUT 
scale, as typical for models inspired by the minimal supergravity (mSUGRA) scenario, RGE 
running will lead to growing flavour- violating effects at low energies. These off-diagonal 
elements can be estimated to have at most the size of [3§1 EE] 

~ " ta (7) (22) 

where ao is an C(l) constant, yo ~ 10~ 5 is the largest Dirac Yukawa coupling (cf. Sec. 14. 7p . 
ms is the universal scalar mass (taken to be some typical superparticle mass), and U rs are 
elements of the PMNS matrix. Then, the maximal (if not even overestimated) value for 
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the branching ratio of [i — > cy can only be 



) 2 

G 2 F m% ~ (m 5 [GeV]) 4 ' 



a 3 (Am?)' in-16 
BrG* e 7 ) ^ < (23) 



where we have used UikUjk ~ 1 and ln(/i //i) ~ 10. However, even for scalar masses of 
only a few 100 GeV, the branching ratio is well below the current 90% CL. limit of the 
MEGA experiment, Br(/i + -> e + 7) < 1.2 ■ 10~ n [571 EHjJH Other channels for LFV are less 
constrained [60], and they are not expected to yield more stringent limits. 
In conclusion, LFV processes make no problems in our case, since the Dirac Yukawa 
coupling yD in the models under consideration is very small, which we could also have 
concluded directly from the fact that we can easily ignore any running, cf. Sec. 14.71 As 
long as one does not include yet another potentially problematic high energy sector, our 
models are safe from this side. 



4.9 Proton decay in SU(5) GUTs 

There is one final remark that is important for GUT theories: In most extensions of the 
Standard Model strong constraints arise from the requirements of gauge coupling unifica- 
tion and of perturbativity, like it was discussed in Ref . [6T] • In particular, the scale M G ut 
at which the gauge coupling unification occurs is related to the proton life-time r p through 
dimension-six operators. In supersymmetric theories, however, the leading contributions to 
proton decay come from dimension-five operators resulting from the exchange of coloured 
higgsinos and winos. These diagrams produce an extremely fast proton decay and they 
have been used to constraint SUSY GUT models. 

The supersymmetric SU(5) in its minimal version has been tightly constrained (if not 
even ruled out) in Ref. [62] assuming that the gauge coupling unification is satisfied and 
imposing the limits provided by the Super-Kamiokande detector, which are particularly 
strong for the p — >■ K + u (t p > 6.7 x 10 32 years at 90% C.L. [63]) and p — > e + 7i° channels 
(r p > 8.2 x 10 33 years at 90% C.L. [63]). 

It is important, however, to note that even if the minimal SUSY SU(5) is highly dis- 
favoured, this does not mean that all possible non-minimal models are excluded as well, 
and some of them have ways to cure the proton decay problem: There exist several works 
in which non-minimal SUSY SU(5) models have been vastly analyzed with the purpose 
of solving the proton decay problem. A possible way out is achieved by a more elaborate 
Higgs sector such that the mass of the Higgs triplet can be pushed to very heavy values, 
suppressing in this way dimension five operators. This is the so-called doublet-triplet split- 
ting problem, see, e.g., Ref. [65] for a review. In this framework, however, it can be difficult 
to suppress proton decay and at the same time not to spoil gauge coupling unification. 
Other ways to suppress or eliminate completely dimension-five operators have been ana- 
lyzed in the context of extra-dimension models, see Refs. [66] [67] 165] for five-dimensional 

3 Note that the next generation experiment MEG currently provides a slightly worse value, Br(/i + — > 
e+ l) < 2.8 • 10~ n at 90% C.L. [59], which is expected to improve considerably within the next few years. 
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SUSY SU (5) theories, and in the context of a flipped SU(5) model, where the up and down 
Yukawa couplings are reversed with respect to the standard 577(5), see Ref. [69J. Flavour 
symmetries have also been used to suppress dangerous proton decay operators [70] . More- 
over, we want to remind the reader that the proton can naturally become almost stable if 
Mgut does not concide with the gauge coupling unification scale [71] . Finally, small group 
theoretical factors can help |72j . 



5 SU(5) inspired models with two FN fields 

We now want to give explicit working examples that yield a simultaneous explanation of 
the low-energy neutrino and charged lepton data, as well as a working scenario with one 
keV-mass sterile neutrino. To do so, we start with an extension of the model presented 
in Ref. [26], which was based on Ref. [13] and yields reasonable agreement with data. We 
will then relax three of their assumptions, in order to be able to construct fully working 
models: We will modify the charge assignments for the right-handed neutrinos, the ones 
for the left-handed lepton doublets, and we will depart from the fully democratic structure 
of the Yukawa matrices. Furthermore, we will present explicit analytical and numerical 
results in Sees. 15.11 and 15.21 

Let us start by recalling the ingredients taken from Ref. [26]: First of all, the model needs 
two FN flavon fields 01,2 which obtain complex VEVs. Physically, we can always choose 
one VEV to be real, in our case Then, the decisive quantities are the ratio A between 

the VEV (0i) and the high-energy scale A, as well as the (complex) ratio R between the 
VEVs: 

A=M, R = <& = Ro e»«, (24) 

where i?o and ao are real numbers. Furthermore, as explained in Ref. [13], it is also 
necessary to introduce an auxiliary Z 2 symmetry, in order for the phase a$ in Eq. (I24p to 
finally be responsible for CP violation. Next, we adopt the following FN charge and Z 2 
assignments inspired by Refs. [H] and [13], respectively, for the FN fields and the lepton 
doublets, as well as for the right-handed charged leptons and neutrinos: 



e lj2 

£l,2,3 
#1,2,3 



" 1 1 5 



(-1,-1; 

(a + 1, a, a; +, +, — ), 
(3,2,0;+,+,-), 

(91,92,93;+,+,-), (25) 



where a = 0, 1 (see Tab. H]). In general, we denote the FN charges of the flavon fields by 
(61,62), the ones of the lepton doublets by and the ones of the right-handed charged 
leptons by hi, as in Sec. |2j The charges gi of the right-handed neutrinos will be chosen 
according to our Scenarios A and B introduced in Sec. EJ The key point is that the FN 
charges of the right-handed neutrinos will drop out of the light neutrino mass matrix when 
constructing the light neutrino mass matrix using the seesaw mechanism of type I or II, 
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Field 


In 


u 


u 


Assignment 1 (a = 0) 


(1,+) 


(0,+) 


(o,-) 


Assignment 2 (a = 1) 


(2,+) 


(1,+) 


(i,-) 



Table 4: Family and Z 2 charges for the lepton doublets. 



as explained in Refs. j26j HTJ |43]. This allows us to freely choose gt without changing the 
appearance of the light neutrino mass matrix. It may be, however, that the charges gi are 
constrained by other consistency conditions (cf. Sec. 14.5ft . 

Before proceeding, we also want to show how the mass matrices are constructed. The 
most general (seesaw type II) Lagrangian that leads to masses in the lepton sector is 



a+b=ki+fj 

C = - E Y2 j e^HL jL X^ + h.c. • 

a,b,i,j 
a+b=fi+fj 

~ E ^Lf<L jL X a X + h.c. 

a,b,i,j 

where, by using Eq. (124"1) . 



a+b=gi+fj 

Y^N^HL jL X a X + h.c. (26) 

a,b,i,j 
a+b=gi+gj 

E if- R M%{N jR ) c XlX\ + h.c., 

a,b,i,j 



A ^"(t)° (t)' , = xa+bRb - ( 2? ) 

The matrices Y e , Yd, and M# are the charged lepton Yukawa matrix, the Dirac neutrino 
Yukawa matrix, and the uncorrected right-handed Majorana neutrino mass matrix, respec- 
tively, which where all taken in Ref. [26] to have a democratic structure. In addition to 
that, ffii is the uncorrected left-handed Majorana neutrino mass matrix, which is present 
in type II seesaw cases. Note that the matrix elements of Mr and rhi are all of the same 
order, i.e., not yet corrected by FN contributions. The sums run over all possible values 
of a and b that fulfill the condition of full cancellation of the FN charges in each term. 
Furthermore, certain terms may violate the Z 2 parity and must hence be set to zero. 
Using all this, we can derive the mass matrices for four different cases: Assignment 1 or 2, 
each combined with Scenario A or B (for type I or type II seesaw), respectively, always 
separated by commas in the respective equations. For the charged leptons, we obtain 

/Y f } 1 B 2 4A 3 ' 4 Y e l2 B 2 X 2 > 3 Y e 13 B 02 RX 2 > 3 \ 
M^ = v \Y e 21 B 2 X 2 ' 3 Y 22 B 2A 1 ' 2 Y 23 RX^ 2 . (28) 
VYftfA 1 ' 2 0, Y e 32 RX YfA - 1 / 

The right-handed neutrino mass matrices for Scenarios A and B turn out to be 

(M R l B && X^ M R 2 B 24 X 3 ' 5 M R 3 RB 2 X 3 > 4 \ 
M^' B) = M 22 B , 2 X ' 2 0, M 23 RX , (29) 

V Mf J 
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while the left-handed ones for Assignments 1 and 2 are given by 

^mi 1 5 24 A 2 - 4 mfB 02 X 1 ' 3 mfRB 02 \ 1 ^ 
m^ 2) = I •' mfBo >2 \°> 2 0, mfRX 2 (30) 



™33r> \0,2 



Finally, the Dirac mass matrices for Assignment 1 are 



'YgB^* Yb 2 B 2A \*<* Y™RB 2 \ 



3,4\ 



m D 



(1A,1B) _ „, | ^1^1,2 y^2 A 0,l Qj y ^ 3jRA | ^ (31) 



Yj?R\ y33 



while the ones for Assignment 2 turn out to be 



m D 



YgB i6 \ 5 >« YgBJ?* Y^ 3 RB 2 4 A 4 ' 5N 

(2A,2B) = ^ [ ygX^a y^ff A 1 ' 2 Y^A 1 '* I , (32) 



y^EA 2 yj 2 i?A yj 3 a 



where u = 174 GeV is the electroweak VEV, and B 2n = 1 + R 2 + ... + i? 2n . Note that, due 
to #2,3 = in Scenario A, the 23-entry of Mr is actually forbidden by the Z 2 parity in that 
case. 

From Mr in Eq. ( l29l . one can immediately calculate the mass eigenvalues for the right- 
handed neutrinos as functions of the right-handed mass scale Mq\ 



A(3,0,0): 



B(4,l,0): 



Mi 


= M A 6 


2^1- 


f Rl + 2# 2 cos(2a ) , 


M 2 


= M , 






M 3 


= M (1 


+ A 6 [14 


- i?o(3 cos(2ct ) + 3i?Q cos(4a ) + i? 4 cos(6a )]) 


Mi 


= M A 8 


2^1- 


f Rl - 2i? 4 cos(4a ) , 


M 2 


= M A 2 






M 3 


= M ( 


1 + i? 2 A 2 


cos(2a )) . 



Indeed, the eigenvalues of Mr show just the hierarchical structure that we have expected: 
If we suppose that Mi is of O(keV), then M should be about 10 6 keV ~ 1 GeV for 
Scenario A, or about 10 8 keV ~ 100 GeV for Scenario B. In any case, we would have a 
low-scale seesaw to work. One can of course raise the possible value for M by simply 
assigning an even higher charge g\ to the first generation right-handed neutrino: Already 
for <7i = 5, one could have M ~ 10 TeV, while one would only have to be careful to keep 
9i — 92 > 3, as explained in Sec. |3j 



5.1 Analytical results 

In this section, we will exemplify for Assignment 1 (and seesaw type I), how one can 
arrive at analytical approximations for the masses and for the PMNS-matrix. We will also 
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see explicitly that it is necessary to depart from the democratic structure of the Yukawa 
matrices in order to obtain sensible results (cf. Sec. 14 . 6 j) . i.e., the FN mechanism alone is 
not strong enough to fully explain the data. It will turn out, however, that this simple 
analytical consideration is still not perfect, even in the non-democratic case, and we have to 
rely on numerics in order to find quasi-perfect models, which will be discussed in Sec. 15.21 



5.1.1 Democratic matrices 



The natural starting point is having democratic forms of all Yukawa and bare mass matri- 
ces, i.e., 

Y? = y e , Yg = y D , and M% = M Vi, j. (33) 

In this case, the type I light neutrino mass matrices can be easily calculated from Eqs. (I2"9"j) 
and (I3ip . using the seesaw formula m u = —m^M^m^. For Scenario A, this results in 



m 



(1A),I 



/A 2 (-2B 2 B4(R 2 +l)+Sl+S 6 R 2 +S 6 ) 

V 



B'j(R?+l)-B e 

-X 
-RX 



whereas the corresponding expression for Scenario B is given by 



-A 


-RX^ 


-1 








-1 / 


by 





(34) 



(1B),I 



m u0 



(-B 2 X 2 -X 

(-BlR 2 +Bl(B2-R 2 -2)+2B 2 B A R 2 +B fi ) 

\ -RX 



B%R 2 ~2B 2 B4R 2 -B 2 B 8 +B 2 +B & R 2 





-RX\ 


-1 / 



(35) 



where m u0 = ^j^- in both cases. In order to obtain the corresponding mixing, we need to 
diagonalize the two neutrino mass matrices in Eqs. ( 1341) and (13"5"|) . as well as the charged 
lepton mass matrix for Assignment 1 from Eq. (I28p . It turns out that all these matrices can 



be approximately diagonalized by applying a series of small and large rotation matrices, all 
of which are (at least approximately) unitary. Such a stepwise diagonalization, although 
certainly not suitable for general mass matrices, is ideally suited for FN models, since the 
corresponding matrices intrinsically involve the small parameter A in which all rotation 
matrices can be Taylor expanded. Note that, in order to derive the leptonic mixing matri- 
ces, we follow the conventions used in the Mixing Parameter Tools (MPT) package [50], [73]. 
which will also be used later on for our numerical computations. 

In order to obtain the (unitary) charged lepton mixing matrix U e , we have to diagonalize 
the squared charged lepton mass matrix M\M e by 



In our case, Eq. 



UlM\M e XJ e = diag(m2, m% m 2 T ). (36) 
, it turns out that the charged lepton mixing matrix is given by 

U e = U A U B U c U D U E U F U Gl (37) 
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where the respective matrices are all at least approximately unitary and are reported in 
Eq. (lA-ip . These subsequent transformations bring M\M e to an approximately diagonal 
form, from which one can read off the expressions for the charged lepton masses to be 

m e = m A 3 Rq , 



m T 



m A (l + A 2 cos(2a ) + fl °~f +3 ] ) 
m (1 + §i? 2 A 2 ) . 



From these relations, we can determine the mass ratios: mg/m^ ~ -RqA 2 and m^/m,- ~ A. 
Using the measured values of m e , m M , and m r , we find the sizes of the parameters to be 
roughly 

A ~ 0.06 and R ~ 1.18 . (38) 

As noted in Ref. [39] , the most advantageous choice of the parameter A turns out to be a bit 
below the standard choice 0.22. The phase a u is not fixed by the lowest order expressions 
in Eq. ( |38|) . but one can choose it to have the value «o = 0.67 in order to make the 
0(A 3 )-correction to m At /m r vanish]^ 

Let us now have a look at the neutrino mass matrix for Scenario A. The matrix m^"^' 1 in 
Eq. (1M|) can be diagonalized by a unitary matrix U v = uj) 1 ^' 1 according to U^m^ U v = 
diag(mi, m,2, m^), where the mass eigenvalues can still contain complex phases. Also 
this can be done by a stepwise diagonalization resulting in 

U v = U a U^U s U e U c U n , (39) 

with the respective pieces being all at least approximately unitary and are given by 
Eq. (O). Note that the purpose of the last matrix U v is only to reshuffle the mass 
eigenvalues in order to accommodate for the correct ordering, since the resulting pattern 
(two larger eigenvalues and one smaller one) can only be realized in inverted ordering. It 
is exactly this point, which will change in the non-democratic cases, and this is also the 
reason why the democratic cases provide a worse match to the data. The neutrino mass 
eigenvalues can be determined as the absolute values of the diagonal entries of the resulting 
mass matrix. They turn out to be 

\m x | = m u0 [1 + A 2 (l + Rl - y /(-l + R 2 ) 2 + R 2 cos 2 q )] + G(A 4 ) , 
|m 2 | = m„ [1 + A 2 (l + Rl + y/(-l + Rl) 2 + R 2 cos 2 a )] + 0(X 4 ) , 
|m 3 | = m^pA 2 . „ $ +0(A 4 ). 

1 ' 2^+2^ cos(2a )+l V ' 

The mass square differences are, at lowest order, 

Am 2 Q - 4m 2 A 2 y/(-l + R 2 ) 2 + R 2 cos 2 ^ + 0(A 4 ) , 
Ami - m 2 (1 + 2A 2 [1 + R 2 - y/(-l + R 2 ) 2 + R 2 cos 2 a ] + 0(X 4 ) . 



4 For Assignment 2, we would have obtained m e /m fl ~ Rq \J\ + RqX 2 and m ll /m T ~ , leading 
to A ~ 0.10 and Rq — 0.80. The phase «o is unconstrained, but the choice ao = 0.60 turns out to be 
numerically convenient. 
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Using the numbers from Eq. fl38l) . one can predict the ratio =^ - [0-018,0.005], respec- 



tively for a = [0, 7r/2], where the choice «o = 0.67 leads to the value 0.015. To obtain a 
fair agreement with the measured value of about 0.031 we would need cto = and thus no 
CP violation. In addition to the imperfect prediction for , the democratic form of the 
Yukawa matrices does lead to a bad mixing. This fact is visible when looking at the full 
PMNS-matrix, C/pmns = U^,U U , which is given by 



/ o 



u 



(1A),I 
PMNS 







-Ci 



Rq 



(1A),I A 2 



Ci 



u 



31 



Ro 



+ u 



(1A),I A 2 



U 



21 

(1A),I^2 
32 



l + (i$ + l)A 2 \ 




A 2 



(9(A 3 ), (40) 



/ 



where 



U (1A),I 
"21 



-Rq(2 cos«o + i sin a — 2£i] 



it 



2^ + Ci 2 



(1A),I 
31 



Ro (2-Rq + (i [i sin a — 2 cos «o] ) 



u 



(ia),i _ (2#o + Ci [2 cos a - i sin a ] ) 



32 



2^ + Ci 2 



(41) 



and (j is defined in Eq. (1A-4I) . Obviously, this matrix does not have the desired form, since 
the angle # 13 is extremely close to the maximal value 7r/2, instead of being small or even 
vanishing. 

Let us now check if Scenario B can give a better match with data. Here, the unitary 
matrix U u = Uv 1 is given by 



U v = U'JJ'M'U'sim , 



(42) 



where the individual rotation matrices are reported in Eq. ( 1A-5I) . Again, the last matrix 
corrects for inverted ordering. The eigenvalues read 



I will 
\m 2 \ 

\m 3 \ 



m uQ (1 + 2i? 2 A 2 ) + 0(X 4 ) , 

(R*(\ 2 +2)-2R2\ 2 cos(2a )+A 2 ) 

wi^o 



m u0 , 



Rl^R^R 2 cos(2a (1 )+l 
^ 2 v / J Rg+2flgcos(2«o)+l + 



C(A 4 ) 



and the mass square differences are, at lowest order, 



Am 2 A 



lA.) 



4Rg 



m 



v0 



(-4 (Ro - 1) A 2 + ^ 
0(A 4 ) . 



4i? 2 A^ 



Then, even with the most suitable value of a = ~ , the ratio of mass square differences, 
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^ m jP ~ 0.35, turns out to be much too large. Also the PMNS matrix, 



/ 



(1B),I 
PMNS 



U 





(lB),Ix 2 
21 A 



K R l- 



-A 



1 + 



i? 



A 2 



^0 



V (1 + ^oA 2 ) 



2ie' 



A 



where 



(1B),I 



3i?^+2cos(2a )-Ro- 1 

3Rq 







-C(A 3 ) 

(43) 
(44) 



2cos(2a ) J R 2 l - 1 

turns out not to be much better than the one from Eq. ( 140~1) . From the above calculations it 
is clear that the hypothesis of democratic Yukawa matrices does not lead to a satisfactory 
PMNS matrix. In the following section, we will analyze the non-democratic case to see if 
we might be able to achieve better agreement with the neutrino data. 



5.1.2 Slightly non-democratic matrices 

The problem of fully democratic matrices not being perfectly suitable for 577(5) inspired 
FN models has already been discussed in Ref. [H]. The key point to obtain a better 
agreement with data is to change the form of the light neutrino mass matrix. Ideally, they 
should have the form 

'A 2 A A\ 

| A 5 2 , (45) 
A 1/ 



with 5 2 
form 



Ami 
Am 2 A 



0.18 > A, whereas our matrices in Eqs. ( |34"|) and (135]) instead have the 



(46) 




which looks similar, but causes a major difference for the mass eigenvalues and for the 
mixings. A suitable modification can be obtained by simply choosing the Yukawa couplings 
Yp 2 = Y^ 2 = 5 y D instead of Y^ 2 = Y£ 2 = y D . Indeed, due to 5 ~ 0.30, this is just a 
relatively mild deviation from the fully democratic case, and it can easily be justified 
by varying the entries in the Yukawa matrices in a certain range around their natural 
values %)£>. 

Using this modification, the new light neutrino masses look like 



/ \ 2 (-2B 2 B 4 (R 2 +l)+Bj+B e R 2 +B 6 



III 



(1A),I 



non— d. 



WW 



i)- 
-R\ 



B 6 



-s x 

-S 2 





"I / 



(47) 



22 



m (1B),I l h 

" b v non— d. 



-B 2 X 2 
\ -RX 



-b 2 r 2 +b 2 (b 2 



S \ 

-R 2 -2)+2B 2 B4,R 2 - 



B%R 2 -2B 2 B i R 2 -B 2 B a +B 2 +B 8 R 2 





-i?A\ 


"I / 



(48) 



The matrix in Eq. ( 14"T|) can again be diagonalized by a unitary matrix t/^, which is now 
given by 

U v = U a UpUJJsU t , (49) 

with the respective pieces reported in Eq. ( ]A-7j) . Note that now, there is no need to correct 
the mass ordering, since the ^-term lowers the second eigenvalue down such that the mass 
square difference between the first and the second squared eigenvalues can now be as small 
as the measured value of AttIq. The neutrino mass eigenvalues can again be determined 
as the absolute values of the diagonal entries of the resulting mass matrix. They turn out 
to be 



II \2 H n 

mi = m uQ A — ; . „° 

11 2 v /R*+2/?2 C os(2«o)+l 

|m 2 | = m uQ (5 2 + 2A 2 ) + 0(A 4 ) , 

|m 3 | = m uQ (1 + 2R 2 X 2 ) + 0{X A ) 



0(X 4 



The mass square differences are, at lowest order, 



Am.Q 
Am 2 A 



™l 5 2 (5 2 + 4A 2 ) + 0(X 4 
m 2 (1 + ARIX 2 ) + G(A 4 



Using the numbers from Eq. ( 138]) . as well ad 5 2 ~ 0.18, we now predict the ratio 



Am? 



0.034, which is a quasi perfect match to the data. In fact, the trick is that the ratio between 
the mass square differences is just given by <5 4 , which precisely justifies the choice made for 



5q. The PMNS-matrix turns out to be 



U, 



(1A),I 



PMNS 



(\ 



\ 



(Rl + 1) * 

(<5 -l)A 
<5o 





iRoe 



i(6 -l)\ 
So 

-l+S +e 2la 





(1A),I X 2 
23 A 



A 2 z(l + 



(9(A 3 ), (50) 



R 2 X 2 ) j 



where u 



(1A),I 



Rn 



i(5%+5o— 1) cosa 
-1 



Although this matrix is still not 



*23 - ^0 ^ ^Zi I" Sin «0 - 5 °^l"° 

really perfect, the value of the angle 9 is is much better than before, since 13 ~ 0(X 3 ). 
Hence, we can be optimistic to find quasi-perfect numerical models when perturbing the 
non-democratic mass matrices. 

Finally, for the non-democratic Scenario B, we have a neutrino mixing matrix given by 



U' a U'^U' s U' £ , 



(51) 



23 



with the respective pieces given in Eq. f ]A-9j) . Also in this case, there is no need to correct 
for inverted ordering. The eigenvalues read 



| mi | 

\m 2 \ 
\mz\ 



m u0 , 



2 v / J Rg+2i?2cos(2Q U )+l rn(x ±\ 

X 2F% + U \ X ) 

[R±(25 2 +\ 2 )-2R 2 \ 2 cos(2a )+A 2 

RI^/rJ^rJ cos(2q ) + 1 



m u0 (1 + 2R 2 \ 2 ) + C(A 4 ) 



+ 0(X 4 



and the mass square differences are, at lowest order, 



Am| ~ 4m 2 5 2 ^2 



Am 2 A 



mto (1 



AR 2 \ 2 



R%-2R%cos(2a )+l 
0(A 4 ) . 



+ o(\ 4 



' Am 2 A ' 



Again taking the most suitable value of ao = § > ^ ne ra ti° °f mass square differences 
is now 0.046, which is still too large but much closer to the actual value than before. The 
PMNS matrix is given by 



/1 



U 



(1B),I 
PMNS 



where 



1 

2<5 



2<5 ^ 







2 i((28 -l)R 2 +e 2 ^o) 
A 27^ A 



II 



[Rl-e 2ia o 

2Rp 
(1B),I A 2 



32 



A 5 





23 A 



+ 0(A 3 



z(l + ^A 2 )7 



(52) 



(1B),I 



23 



n 2 3 

(1B),I 
32 



™32 



i? -ie~ ia ° + 



™23 



(4<5 4 -1)^ + 2 cos(2a ) J Ro ~ 1 , 
5 (2e~ 2iao 5%R 2 + R 2 - e 2ia °) ((R 2 + l) sina - i (R 2 - l) cosa ) 

ie- Mao R n 32 
(45 4 - 1) i? 4 + 2 cos(2a ) J Ro 1 1 ' 



-2e 6iQ0 i^ 



i? 2 (2i^ 3 - 8 + 1) + 



+e 2 ^((45 4 + 5 -l) i2j + *o-l)- 



(53) 



Also in this case, the value of the angle #13 is much better than in the democratic case, since 
6* 13 ~ 0(X 3 ). However, to find quasi-perfect models, we will have to turn to a numerical 



study, which is presented in the next section. 



5.2 Numerical analysis 

Finally, we want to present a few quasi-perfect models and their predictions. Usually, and 
in particular in the context of FN inspired models, so-called scatter plots are presented in 
the literature. These plots are supposed to indicate that the model is consistent with data 
in a considerable region of the parameter space. However, we do not consider this approach 
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1AI 



M 



R 



Matrix 



M, 



eO 



0.815 2 A 3 
2.005 2 A 2 
3.71i?A 



l.UB 2 X 2 

1.13A 





0.29RX 2 

2.50RX 

0.35 



0.7554A 4 
0.51A 
3.32i?A 



0.155 2 A 3 

0.13 





0.38£ 6 A b 
0.315 2 A 3 
1.265, R\ 3 



0.31S 2 A 3 

4.18 






1.26B 2 Ry 


4.81 



Table 5: Mass matrices of model 1AI. The prefactors parametrize the mass scales imposed 
on the concrete model at a certain energy scale. 



1BI 



m D 



M 



R 



Matrix 



0.91£ 2 A 3 
3.805 2 A 2 
0.79i?A 



2.265 2 A 2 

2.51A 





4.33.RA 2 

3.63.RA 

0.15 



3.205 4 A 5 
2.275 2 A 2 
1.69RX 



O.I554A 4 

0.030A 





3.73E 2j RA 4 
1.23.RA 
0.66 



O.33.08A 8 
M ( 2.6854A 5 
1.63S9.RA 4 



2.68B4A 5 
2.10S2A 2 
0.83i?A 



1.63E2.RA 4 

0.83i?A 

0.49 



Table 6: Mass matrices of model 1BI. The prefactors parametrize the mass scales imposed 
on the concrete model at a certain energy scale. 



as too useful, since first such results always depend on the statistical measure used in the 
generation of the random numbers involved, and second because there is essentially no 
method to determine if indeed a large part of the parameter space is investigated or rather 
only a tiny patch, due to the gigantic complexity of higher- dimensional non-Cartesian 
spaces. Furthermore, there is no way to decide whether a model can be considered as 
"good" or "bad", just because a certain fraction of random parameter choices leads to 
compatibility with data. 

We will rather take on a contrary approach and try to find four fully working examples by 
perturbing the mass matrices from Eqs. ff28l) . fT29|) . fl3Tl) . and fl32|) . where the Dirac Yukawa 
coupling matrices Yd are taken to be non- democratic, according to Sec. 15.1.21 If these 
examples are in agreement with data, we consider the FN approach as being predictive in 
the sense that small departures from the analytical forms of the mass matrices are perfectly 
enough to be in full agreement with low-energy neutrino data. 

Based on the analytical results from Sec. 15.1-H we choose the parameters A, Rq, and a>o 
in the following way: 
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Matrix 
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M e0 


/ / ±.yO-D4A 0.IV-D2A O.UO£)2-UA \ 

4.755 2 A 3 1.5152A 2 3A2RX 2 
\ 3.17RX 2 0.70RX 0.28A / 




m D 


moo j 


' 3.79B A X 5 0.00575 4 A 4 4.605 2 -RA 4 

0.6852A 2 0.097A 1.20.RA 
^ 1.12RX 2 3.3CLRA 4.85A 


) 


M R 


M | 


' 3.02£ 6 A 6 1.2852A 3 0.9W 2 RX 3 \ 

1.28B 2 X 3 4.32 
v 0.935 2 -RA 3 3.08 J 





Table 7: Mass matrices of model 2AI. The prefactors parametrize the mass scales imposed 
on the concrete model at a certain energy scale. 

• Assignment 1: A = 0.06, R = 1.18, a = 0.67, 

• Assignment 2: A = 0.10, R = 0.80, a = 0.60. 

The next step is to generate a number of perturbed charged lepton matrices according to 
Eq. (128]) . and check which of them yield a quasi perfect prediction (±1%) for the charged 
lepton mass ratios m e /m M and m^/m T . We furthermore generate right-handed neutrino 
mass matrices according to Eq. (129]) , whose smallest eigenvalue should always be at least six 
orders of magnitude smaller than the one of the second to lightest mass eigenstate. We then 
use these matrices together with the Dirac mass matrices from Eqs. ( 13~TT) and ( 1321) to gen- 
erate light neutrino mass matrices according to the seesaw formula, m v = —m^M^mu, 
which yield a quasi-perfect prediction (±1%) for the ratio of neutrino mass square dif- 
ferences, Amg/Am 2 ^. Finally we combine the lists of matrices and calculate the mixing 
parameters and phases using MPT [501 173] . All models that lead to predictions of the 
mixing angles that are consistent with the current 3a ranges [19] will survive the test. 
The corresponding mass matrices of the quasi-perfect numerical models for all four com- 
binations of Scenarios A and B, as well as of Assignments 1 and 2 can be found in Tabs. [5] 
to [HJ Indeed, all mass matrices have just the structure that we have desired: Only slight 
perturbations of the democratic form, with the exception of the Dirac mass matrices that 
have smaller 12 and 22 entries. Accordingly, as we had expected, the FN models presented 
here are predictive up to the exact values of the mixing angles and CP phases. In particu- 
lar, all charged lepton and light neutrino mass matrices lead to the correct ratios of masses 
or mass square differences, respectively. 

The numerical predictions are summarized in Tab. [9], and the mixing angles and masses 
are also displayed in Fig. [3J As can be seen, all mixing angle predictions are in agreement 
with the 3a ranges of the oscillation parameters , because the models have been selected 
in that way. Furthermore, as explained in Sec. 14.71 these predictions are effectively scale 
invariant. However, a non-trivial prediction is the actual values of the masses (in terms of 
the absolute neutrino mass scale m u , which is determined by the other scales involved), and 
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2BI || Matrix 

| 7 2.71/3 4 A 4 4.17t3 2 A 3 3.74/3 2 i?A 3 

M e M e0 1.39/3 2 A 3 0.19/3 2 A 2 3.19.RA 2 

V 0.9LRA 2 0.38#A 0.40A 

7 1.57/3 6 A 6 0.078 /3 4 A 5 1.795 4j RA 5 

?77, m D0 4.87/3 2 A 3 0.1673 2 A 2 4.85i?A 2 

V 3.47i?A 2 LS2RX 1.25A 

/ 0.61/3 8 A 8 1.12/3 4 A 5 0.53/3 2j RA 4 

M r M 1.1273 4 A 5 3.85/3 2 A 2 4.33.RA 

I \ 0.53/3 2 #A 4 4.33/ZA 3.21 

Table 8: Mass matrices of model 2BI. The prefactors parametrize the mass scales imposed 
on the concrete model at a certain energy scale. 
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m z jm v 


Hier. 


1AI 


0.28 


0.018 


0.54 


4.21 


0.58 


1.25 


0.0014 


0.19 


1.03 


NH 


1BI 


0.31 


0.035 


0.41 


4.98 


2.35 


2.89 


0.099 


0.59 


0.0028 


IH 


2AI 


0.31 


0.0015 


0.40 


3.88 


2.27 


0.54 


1.7 • 10" 8 


0.020 


0.11 


NH 


2BI 


0.27 


0.021 


0.57 


3.83 


1.96 


2.71 


0.0018 


0.0054 


0.029 


NH 



Table 9: Predictions of the four numerical models. "Mod." stands for "model", s 2 - = 
sin 2 #jj, and "Hier." stands for "hierarchy". 



in particular the dominant prediction of normal mass ordering (or, due to the hierarchical 
structure of FN matrices, normal hierarchy), which we had already anticipated in Sec. 15.1.21 

6 Conclusions 

In this paper, we have shown how to carefully derive models explaining the appearance of 
one right-handed (sterile) neutrino that has a mass at the keV scale, while at the same time 
predicting leptonic mass ranges and mixing parameters in full agreement with experiments. 
These models were based on the famous Froggatt-Nielsen mechanism, which is a well-known 
possibility to create strong hierarchies between fermion masses. To our knowledge, this 
is the third known type of models that can successfully explain the existence of a keV 
sterile neutrino Dark Matter particle, the first two being based on soft breaking of flavour 
symmetries [HI [T5] and on Randall- Sundrum warping [T7], respectively, where the former 
exploits symmetry breaking effects to lift a massless state to the keV scale while the latter 
uses the exponential suppression of UV-brane physics to strongly suppress the natural 
right-handed neutrino mass scale. We have instead made use of the possibility to assign 
FN charges to the different fermions in such a way as to suppress certain mass matrix 
elements, whose sizes are reduced by roughly one order of magnitude per unit FN charge. 
In this way, it is possible to suppress one right-handed neutrino mass strongly enough to 
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Figure 3: The predictions of the numerical models for the leptonic mixing angles as for the 
masses, the latter in units of the light neutrino mass scale m v . 



be around the keV scale, while the other two can easily be considerably heavier (a mass of 
about GeV is the lower bound to be fulfilled, but they could easily be much heavier). One 
important bonus of this approach is that, due to the structure of the FN charge assignments, 
the seesaw mechanism will be guaranteed to work if it works with all FN charges set to zero, 
so that we can be sure that the existence of a keV particle will not lead to any problems 
from that side. Furthermore, we discuss the requirements and potential problems that could 
arise when applying the FN mechanism in certain frameworks. Interestingly, although FN 
charge assignments might seem relatively arbitrary at the first sight, it is easy to find 
situations where the different conditions are so restrictive as to render the different sectors 
incompatible. This discussion leads to a systematic reduction of the arbitrariness involved, 
and it turns out that an ideal framework for our purpose is an £77(5) GUT, augmented by 
two FN flavon fields. Due to the relatively low seesaw scale involved, typical constraints 
from RGE running and the corresponding LFV effects can be avoided. Finally, we show 
analytically that democratic structures of the mass matrices are not enough to explain all 
mass ratios, but once we depart from this structure only slightly, the agreement with the 
data improves. In order to find fully working models, which are also in agreement with 
the experimental constraints on the leptonic mixing angles, we finally perform a numerical 
analysis that perfectly justifies our argumentations given before. 
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Appendix A: Stepwise diagonalization 

We report in this appendix the explicit expressions of the approximately unitary mixing 
matrices that diagonalize the charged leptons mass matrix and the neutrino mass matrices 
for Assignment 1, Scenarios A and B, with both democratic and non- democratic Dirac 
Yukawa matrices. 



Charged leptons 

Here we report the mixing matrices that diagonalize the charged leptons mass matrix, as 
reported in Eq. ( 1571 ): 



U A 



U 



D 




U G = -gtX 3 1 , (A-l) 



with 



-iot 

> 



a\ = b\ = c\ = Roe' 

e x = Roe* 00 (R 2 e- 2iao - 2R 2 + e~ 2ia ° + 2) , 

h = -R 3 e- ia ° : 

9l = -R 2 (R 2 + e- 2iao ). (A-2) 
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Light neutrinos, Assignment 1, Scenario A, democratic 

Here we report the mixing matrices that diagonalize the light neutrino mass matrix, see 
Eq. (1391) . for Assignment 1, Scenario A, and democratic Yukawa matrices: 
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— i?o + 1 + -Ro) 2 + -^o cos2 sec a . 
Light neutrinos, Assignment 1, Scenario B, democratic 



(A-4) 



Here we report the mixing matrices that diagonalize the light neutrino mass matrix, see 
Eq. (142p . for Assignment 1, Scenario B, and democratic Yukawa matrices: 




u: 
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with 



a[ = Roe-** 
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2Rl 

R cosa [l + 2Rq + - 6R 2 cos(2a )] 
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61 = 8^ • (A " 6) 

Light neutrino, Assignment 1, Scenario A, non-democratic 

Here we report the mixing matrices that diagonalize the light neutrino mass matrix, see 
Eq. (H9l) . for Assignment 1, Scenario A, and non-democratic Yukawa matrices: 
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Light neutrinos, Assignment 1, Scenario B, non-democratic 

Here we report the mixing matrices that diagonalize the light neutrino mass matrix, see 
Eq. (1511) . for Assignment 1, Scenario B, and non-democratic Yukawa matrices: 





K=\ o i o | , u' = ( -(#)*a i o | ,[/;=( z z 7 ;a 2 

-W)*A 3 I , (A- 
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